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NOMENCLATURE 


a  - 

a  - 
n 

A(C,t) 

Ai(5,t) 

b  - 

b  - 
n 

C  - 
P 

Cg  - 

C  - 

c  - 
n 

Cj.Cz  - 

c  - 

cn  - 
cd  - 
d  - 
dn  - 
D(e,T) 
Dn(X)  - 
e  - 

E(5,t) 
E,E(p,r 
F(C,t) 
F(p,r) 
f (^,z) 


wave  amplitude 

see  eq.  (3.1.2) 

see  eq.  (2. 1.8) 

see  eq.  (2.1.52) 

see  eq.  (B.3d) 

see  eq.  (3. 1.3) 

phase  velocity  -  see  eq 


(2.1.36) 


group  velocity  -  see  eq.  (2.1.28) 
see  eq.  (2. 3a. 5) 
see  eq.  (3. 1.4) 
see  eqs.  (B.l),  (B.2) 
see  eq.  (B3.a) 

Jacobian  Ellyptic  Function  -  see  eqs.  (C. 15) 
Jacobian  Ellyptic  Function  -  see  efter  eq.  (C. 12) 
see  eq.  (B.3b) 

Jacobian  Ellyptic  Function  -  see  eqs.  (C.15) 

-  see  eq.  (2. 1. 19) 
see  eq.  (3.1a.l) 
see  eq.  (B3.c) 
see  eq.  (2. 1.23) 

)-  Ellyptic  Integrals  -  see  Byrd  and  Friedman  (1971) 
see  eq.  (2. 1.22) 

-  Ellyptic  Integral  -  see  Byrd  and  Friedman (19 71) 

-  see  eq.  (D.  1) 


gravity  constant 
see  eq.  (2.1.20) 
see  eq.  (2. 3b. 4) 
still  water  depth 
water  depth 
complex  tmlty 
see  eq.  (3.1a.6) 

see  eqs.  (3. 1. 12)  ,  (3. 1. 13)  ,  (3. 1. 14) 
carrier  wave  number  -  see  eq.  (2.1.8) 
wave  number  -  see  eq.  (2. 3a) 
see  after  eq.  (C.18) 
wave  length,  see  eq.  (3.1d.l) 
group  length,  see  eq.  (3. Id. 2) 

'supergroup  length'  -  see  eq.  (3. Id. 3) 

see  eq.  (3.1.21) 

see  eq.  (3.1.17) 

see  eq,  (C.19) 

see  eq.  (C.17) 

see  eq.  (2.2.5) 

see  eq.  (A. 2) 

see  eq.  (2.2.9) 

see  eq.  (B.4) 

see  eqs.  (C.15) 

Jacobian  Ellyptic  Function  -  see  eqs  (C. 14) 


dimensionless  form  of  t. 


Y  -  see  after  eq.  (3.1.1) 

^1*^2' (A..  8),  (A.  9) 

6  -  see  eq.  (C.24) 

e  -  see  after  eq.  (2.1.4) 

4  -  free  water  surface 

4q,5i,C27.?®®  (2.1.7) 

4^j(C,t),  mil,  Jim  -  see  eq.  (2.1.10) 

^20  -  see  eq.  (2. 1. 17) 

6  -  see  eq.  (2. 1.8) 

Xj ,  j=l,9  -  see  eqs.  (2.1.4^^  to  (2.1.51) 
u(C)  -  see  eq.  (3.1.9) 

Pj,  j=l,6  -  see  eqs.  (2.1.43)  to  (2.1.51) 

C  -  see  eq.  (2.1. 15) 

a  -  see  eq.  (2.1.26) 

o  -  see  after  eq.  (2. 3a. 4) 

T  -  see  eq.  (2.1. 15) 

-  see  eq.  (2.  3c. 6) 

-  see  eq.  (C.30) 

<J>  -  Induced  velocity  potential 

‘^o*'^l’*^2  ^9'  (2.1.16) 

(5,t,z)  ,  m)^l,  j>,m  -  see  eq.  (2.1.9) 

'*’l0’**’20  “  ^2.1.11) 

<t>  -  dimensionless  form  of  ())^q  -  see  eq.  (3.2.5) 

'I'  - 


wave  envelope 


dimensionless  form  of  iJj 

in  section  (3.1)  see  eq.  (3.1.6) 

in  section  (3.2)  see  eq.  (3.2.5) 

stream  function 

dimensionless  form  of 

in  section  (3.1)  see  after  eq.  (3.1e.l) 

in  section  (3.2)  see  eq.  (3.2.16) 

carrier  wave  frequency  -  see  eq.  (2.1.8) 

'.rave  frequency  -  see  after  eq.  (2. 3a. 4) 
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1.  INTRODUCTION  AND  ACHIEVEMENTS 


The  shoaling  of  weakly  nonlinear  surface  wave  groups  is  important  to 
the  understanding  of  coastal  wave  climate  and  coastal  flow  regime. 

In  the  past,  most  efforts  concentrated  on  the  equally  Important 
though  simpler  problem  of  shoaling  of  wave-trains  (i.e.  mono chr omat ic 
wave  groups),  for  details  see  Stiassnie  &  Peregrine  (1980), 

The  first  mathematical  formulation  for  shoaling  of  wave-groups 
was  given  by  Djordjevlc'  and  Redekopp  (1978).  This  formulation  is 
limited  to  cases  where  the  water  depth  is  small  compared  to  the  group- 
length.  Equations  suitable  for  water  depths  of  the  order  of  the  group- 
length  are  deduced  in  Peregrine  (1983) ;  combining  the  constant  depth 
model  by  Davey  and  Stewartson  (1974)  and  the  higher-order  model  for 
infinitely  deep  water  by  Dysthe  (1979). 

The  only  available  solutions  are  those  for  the  shoaling  of 
isolated  wave-packets  (solitons) ,  which  were  originally  given  by 
Djordjevlc'  and  Redekopp  in  their  1978  paper.  They  predict  that  a 
sollton  envelope  can  undergo  fission  only  if  it  propagates  into 
deeper  water.  3y  heuristic  assumptions  for  the  evolution  along  the 
slope,  they  also  estimate  the  number  of  solitons  emitted  after  a 
single  sollton  descends  from  a  shallower  shelf.  A  more  recent  study, 
Turpin,  Benmoussa  and  Mei  (1983)  confirms  these  results  qualitatively, 
but  not  quantitatively. 

To  the  best  of  our  knowledge,  no  results  for  shoaling  of  wave 
gr'jups  (i.e.  modulated  wave-trains)  have  been  presented  so  far.  These 
modulated  wave-trains  are  of  particular  importance  since  almost  every 
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wave-traia  will  eventually  become  modulated  due  to  its  intrinsic 
Benjamin-Feir  instability.  The  main  goal  of  the  present  study  is 
to  throw  light  on  the  evolution  during  the  shoaling  of  a  modulated 
wave-train  and  its  influences  on  the  mean  free  surface  and  the  wave- 
induced  mean  flow. 

The  bulk  of  this  work  is  divided  into  two  parts:  the 
derivation  of  the  appropriate  equations  (in  chapter  2),  and  their 
solutions  (in  chapter  3). 

The  use  of  the  REDUCE  2  algebraic  manipulator  enabled  us  to 
derive  evolution  equations  which  are  accurate  to  fourth  order.  These 
equations  are  valid  for  any  water  depth,  (except  the  extremely  shallow) , 
as  long  as  the  slopes  of  the  bottom  are  sufficiently  mild.  Note 
the  nonlinear  effects  on  shoaling  surface  gravity  waves  in  extremely 
shallow  water  are  discussed  in  a  recent  paper  by  Freilich  and  Guza 
(198A) . 

Two  new  results  are  presented  and  discussed  in  this  report.  The 
first  is  an  approximate  analytical  solution  which  provides  detailed 
information  about  the  physical  processes  involved  in  the  shoaling  of 
wave-groups.  The  second  result  is  the  calculation  of  the  induced  mean 
flow  accompanying  an  evelope-soliton  moving  on  infinitely  deep 
water. 

The  attempts  to  include  the  effects  of  randomness  into  our 
models  have  not  yielded  any  wortwhile  results.  It  seems  that  the  very 
recent  work  by  Longuet-Higgins  (1984)  should  provide  a  good  starting 
point  for  any  future  work  in  this  direction. 


2.1  Fourth  Order  Evolution  Equations 


The  third  order  evolution  equations  for  groups  of  water-waves  moving 
over  an  uneven  bottom  were  first  derived  by  Djordjevic'  and  Redekopp 
(1978).  These  equations  are  valid  only  if  the  water  depth  h  is  of  the 
same  order  of  magnitude  as  the  wave  length  S,^.  In  the  present  section 
we  carry  their  derivation  one  step  further,  to  fourth  order  and 
introducing  a  slight  modification,  we  obtain  a  more  general  set  of 
equations  which  is  valid  when  0(h/ij)  >  1.  All  the  algebraic  mani¬ 
pulations  were  performed  utilizing  the  REDUCE  2  algebraic  manipulator, 
see  Hearn  (1973). 

We  consider  the  evolution  of  a  uni-directional  progressive 
grativy  wave  mi'ving  along  the  x  axis  on  the  free  surface  of  a  homogeneous 
liquid  with  depth  h  =  h(x)  varying  in  the  direction  of  the  propagation. 

The  effect  of  surface  tension  is  neglected,  so  the  analysis  applies 
to  gravity  waves  only.  The  fluid  motion  is  irrotational,  thus  there 
exists  a  velocity  potential  (t>(x,z,t)  satisfying  Laplace's  equation 

(ji  +4’  —  0  (2.1.1) 

XX  zz 

where  z  is  the  vertical  coordinate,  and  z  =  0  is  the  undisturbed  free 
surface. 

The  boundary  condition  on  the  bottom  z  =  -h(x)  is 


(j)  =  -h'  (x)  <p. 


(2.1.2) 


and  the  boundary  condition  on  the  free  surface  z  =  C(x,t)  are 
the  kinematic  condition; 


“  ^t  + 


and  the  pressure  condition 


(2.1.3) 


2g?  +  24.  +  4.5  +  (f.^  -  0  (2.1.4) 

V  X  z 

The  situation  where  the  depth  varies  slowly  in  the  direction  of 
propagation  is  considered,  so  that  properties  characterizing  the  wave  will 
change  slowly  as  well.  A  small  nondlmenslonal  parameter  e  that  measures 
the  slope  of  the  wavy  surface  is  introduced,  and  we  define  the  new 
variables 


X 

T  =  e[/ 


where  C  is  the  group  velocity. 


(2.1.5) 


We 


suppose  that  the  depth  changes  on  the  scale  of 
h  =  h(C)  with  the 


so  that 


The  velocity  potential  and  the  free  surface  displacement  are 
expanded  as  follows: 

<!•  =  <t>^(5,x,z)+{(|)j(5,T,z)e^%2<5*^»*>®^^®+--*+‘'-‘'^  (2.1.6) 

?  =  ?jj(5,T)+{Cj(5,T)e^®+42(5,T)e^^®+...+c.c}  (2.1.7) 

where  0  •»  /  k(C)dx  -  ut  (2.1.8) 


and  c.c  means  complex  conjugate. 


where  *  <^£jCC»T,o),  and  *  denotes  the  complex 


conjugate. 


From  (2.1.14)  and  (2.1.15)  it  follows  that  the  order  of  4^ 
and  <j),_  (o)  are  greater  or  equal  to 


•liin  (o)  <  O(e^);  $  O(e^) 


(2.1.16a,b) 


We  Introduce  the  notation 


'o  =  ^  ^20 


(2.1.17) 


The  next  step  of  the  derivation  is  similar  to  that  given  by 
Djordjevic*  and  Redekopp  with  the  difference  that  we  continue  the 
process  to  fourth  order.  Substitution  of  (2.1.6),  (2.1.8)  and  (2.1.9) 
into  Laplace's  equation  (2.1.1)  gives: 


0(e)e^®! 


’^ll  ’  cosh  k  h 


(2.1.18) 


rf,  .  nz-F  T-k  cosh  k(z4-h) 
*21  cosh  k  h 


[z4-h)  sinh  k(z+h)  -  ho  cosh  k(z+h) 
cosh  k  h 


(2.1.19) 


.  X  cosh  k(z+h) 

*31  ■  co7hk  h~ 


{2kh'At(k.'A  -  |rA„)(z+h)>- 


2  cosh  k  h 


g 


(z+h)cosh  k(z+h)  + - r— r^{cr(hk)  '  A  - 

cosh  k  h 


-A--  A _ -  -^}D  •(z+h)sinh  k(z+h)  (2.1.20) 


+  [.  i- .  k* 


A  -  —A  ICfxfM^]  k(z4-h) 

bC"  ^ttt  2C  ^JLCz+h)  J  ^ 

S  g 


-  ^  /v(z+h)^  ^r**^  + 

C  T '  cosh  k  h 


cosh  k  h 


■*■^2^“  Ukh’)  [3a-2kh(2ci^-l)]  +  kh'a  -  hk‘o- 
g 

-  2kh<J*  — }  A  +  —A  “ 

C  k  T  C  tS 

g  g 


iG. 


±<^  k'  _  T^.  ,  .,,  sinh  k(z+h) 

-  ,  -2  ■; —  D-iDf  -  -r— "-f  (z+h)  - r  r-V — ^ 

k  5  Cg  ^  '  cosh  k  h 

2  C* 

^  k'+2khok'  +  2k  /  IK,- 

8  g 

iha  .  °tt  ik'  ,  ,.,2  cosh  k(z+h) 

”  2^  \Tt*  —  2?-  -  ~  cosh'k  h 


+{  -  ^  (l-kha)A^-ikh'D}  (z+h) 


mlhltj 


^22  '  cosh  2  k  h 


A  .  E(C  T)  -£25lL2kIz±*ll _ I 

^32  cosh  2  k  h  C 

g 


•  (z+h) 


slnh 


2k(z+h) 
cosh  2  k  h 


+ 

(2.1.21) 


(2.1.22) 


(2.1.23) 


where  A,D,G,F,E,  and  M  are  yet  unknown  functions. 


Substitution 

of  eq.  (2.1.6) 

to  (2.1.10)  and  eqs.  (2.1.18)  to 

(2.1.23)  Into  the  expansion  around 

2  =  0  of  the  free  surface  conditions 

(2.1.3)  and  (2.1.4) 

gives : 

and 

r  =  iifl  * 

11  g 

(2. 1.24) 

01^  =  g  k  0 

(2.1.25) 

where 

0  =  tgh(kh) 

(2.1.26) 

^20  ■  g^*^10T“'^20^  g 


(2.1.27) 


_ & 

g 

(2.1.34) 

il-  ^  (l-kh(T)(l-CT-)] 

g 

(2.1.35) 

C 

■  [2  +  (1-ff^)]  ;  C  =  u/k 

g  P 

(2.1.36) 

■  [—  -  12  +  130^-  2a‘»] 

(2.1.37) 

(1-02) 

2 

t 

t 

(2.1.38) 

.  -  .  . ^iuko  ik^  j  *  , 

lojG  +  (1+  ^  )D  +(  „  "  r  ^  ^ 

S’*  S  ' 

(2cr*-l)A  ,+hoiuA,  - 
2CI  ”  ^ 

g 

{3-90^+llcj‘*  -  20®}A^A  + 

2u)0^ 

^  [2kh'  +  k'h  -  20^ (kh)']  A 

(2.1.39) 

'  “  C3kh(3cT‘»+202-5)+9a(l-02)]  A  A^ 

g 

-  (1-02)AD 

(2.1.40) 

=  — [-12kh(l-o2)-404l20]  A  A^ 

4C  03 
g 

(3-0^) AD 

(2.1.41) 

-1 
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>e  ; 


+  U3D^^  +  y4D'^io  ■*'  ^yjlApD  +  UjA^D 


+  yg<|.2oD  +  +  X^A.^^  +  X3|a1^A^+  X^A^A^  + 


+  X5A  +  XgA<t^Q  +  X^A^.^  +  X  A<^  + 

TT  T  C 


+  =0 

Z 

where  yj^,..,,yg  are  given  by  (2.1.32)  to  (2.1.33)  and 

X,  =  -  (2h3k^a’-6h3k3CT5+6h3k3a3-  Zh^k^o  - 

^  401^03 

g 


-  5h^k^(7^+9h2k^CT'*-5h2k2a2+h2k2+  Ahkcr^- 


-4hka3  -  0*^+0^) 


[_3(j+kh(3a2-l)] 
3C- 


(2.1.42) 


2  ..3 


(2.1.43) 

(2.1.44) 


X»  =  -  {7cr2-47a5+48a3-36a+kh  [-7a®+32a^- 

2gC^a3 

-25a‘*-18a2+18]. 


,  ,-ild 

^  2gC„ 


{  -o'^rSo^+khcT  ( l-a^  )  2  } 


X5  =  — (l-kho) 

g 


^6  = 


laikq  21k  ta> 

8^  '  *^8  ’  ec5 


(2.1.45) 


(2.1.46) 


(2.1.47) 


(2.1.48) 


<►'  * 


.  .V 


kv 


Combining  eqs.  (2.1.31)  and  (2.1.40)  and  introducing  the  new  variable 


=  eA  +  e^D  (2.1.52) 

one  finally  obtains  the  modification  to  fourth  order  of  the  cubic 
SchrBdinger  equation 


MfAi  +  1^2^15  ^4^1^10  ^  1*5 1 

TT  “t 


+  ^^1^1  ^  ^31^1!  + 

\  ^TTT  T 

+e"  X^A^^A  +  €X5A  +  X^eA  *  +  eX  A  + 

T  TT  T  T 


+  +  E  XgAj(i>jQ  =  0 

4  z 


(2.1.53) 


On  the  other  hand,  substitution  of  (2.1.18),  (2.1.24)  and  (2.1.11) 
into  (2.1.14)  and  (2.1.15)  and  elimination  of  yields  the  following 


equation  for 


'10’ 


♦10  *  t  no,,-  ^  C  +(1-o=)H|a|I=«>*«A.D)^ 
z  ^  g 

+  2|  ^3  |_  [k|A|2].e3r.  ^  -  |a][AA.^-A.A^j4<)Uk) 

(2.1.54) 


which  can  be  written  to  the  same  order  of  accuracy  as: 

c^4’lo  1  2  2C  r, 


•|^WA.f]-a-iJ^*^-|^KA,A5^^-AtA^,] 


(2.1.55) 
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Laplace's  equation  (2.1.12)  and  the  boundary  conditions  (2.1.2)  at 
z  =  -h  and  eqs.  (2.1.51),  (2.1.53)  at  z  =  0  fom  the  system  of  equations 
from  which  and^j^^  can  be  determined. 

For  the  case  of  infinitely  deep  water,  as  it  will  be  shown  later, 
the  induced  flow  is  of  order  e^.  For  this  case  equations  (2.1.55;  and 
(2.1.53)  become  respectively: 

<^10  (2.1.56) 

z  ®  T 


and 


2aiiC  *^1  1  »  , 

K  .  .  TT  .  2a)k  .  , 

g  1-  g  R&  1  10 


2ik^  .or  T  -^iuk 


The  laist  two  equations  are  identical  to  those  given  in 
K.B.  Dysthe  (1980). 


(2.1.57) 


Note  that  the  following  typographical  errors  were  found  in  Dysthe’ s 
paper: 

His  equations  (2.17)  and  (2.19)  should  be  written  (in  his  notation)  as 
follows; 

r  =  4k‘*[A[2+8ik3(AA*  -  A*A  )-4ik3[A[^  +2(Aik(i  -ii  ) 

3C  X  X  X  Z 

(Dysthe  2.17) 

21  (A  +?iA  )+i5A  -hk  -a|a1^=  -  i(6A  -A  )  + 
t  X  yy  XX  '  ‘  8  xyy  xxx 


+  |i4(AA*-A*A^)-islAl2A^+A(i^-a^) 


(Dysthe  2. 19) 
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2.2  Third  Order  Eyolutton,  Equations 

Equations  (2.1.32)  and  (2.1.55)  constitute  the  following  third 
order  set  of  evolution  equations: 

UjA  +  +  UqA^^  +  V‘4A-<(ijo  +  P51aI^A  +  P^Acfi^o  =  0  (2.2.1) 

^10  2.  2  2C 

^10  +  [-^  +(l-a^)]  ([a[^)  (2.2.2) 

Z  ®  ®  g  T 


Shallow  water 

For  shallow  water,  we  assume  that  the  horizontal  induced  current 

is  uniformly  distributed  with  depth;  and  then  integrating  equation 

(2.1.1)  yields  2a 

-e  fio 

<I>1Q  =  — q7 — ^  (2.2.3) 

z  g 

Substituting  (2.2.3)  for  z=0  in  (2.2.2)  yields 


+(l-a^)](lAl^) 

g  T 


(2. 2. A) 


integrating  (2.2.4)  with  respect  to  x  gives 

♦  -  [-f^ +(1-0^)]  |A|2  +  Q(0  (2.2.5) 

(1  -  1^) 

g 

The  set  (2.2.1)  and  (2.2.5)  is  the  same  as  that  given  by  Djordjevic' 
and  Redekopp  (1978),  except  the  term  ygA({)2Q(?)  which  does  not  appear 
in  DJordJevlc  and  Redekopp. 


9- 


In.  the  sequel  we  vrLll  show  the  necessity  of  including  the  term 
in  the  wave-induced  mean  current  expansion  (see  (2.1.11). 


To  second  order  in  e  the  wave  induced  mean  current  velocity  is  given 


by 


C,(I-  ^ 


(l-a2)]|A|2 

g 


(2.2.6) 


To  find  Q  we  introduce  a  lateral  boundary  condition  of  zero 
averaged  (over  x  )  mass  flow  which  is  appropriate  for  an  impervious 
beach,  as  follows: 

hU+e^k^]A|T  (2.2.7) 


where  the  bar  indicates  the  averaging. 

From  eqs.  (2.p.  77  artd  (2.2.6)  we  obtain 
2u)kC  t#-^ 

Q(0  =  -  (2.2.8) 

(1-  C  7  g 
g 

Integration  of  eq.  (2.2.5)  with  respect  to  t  yields 
1,2  2C 

*  - ^^-r-  [-7^ +(l-a2)]  /  |a|^  dx  +Q(Ot  +  eQ  (O  (2.2.9) 

(1-  -|J)  ^g 

g 

The  first  and  second  terms  in  (2.2.9)  grow  monotinically ,  and 
boundlessley  in  time.  Secular  terras  of  this  nature  are  bound  to  cause 
trouble  in  higher  order  derivation  and  should  be  suppressed.  The  addition  of 
the  term  e;24,2Q(^)t,  to  eq.  (2.1.11)  seems  to  be  the  proper  way  to 
achieve  this  goal.  Substituting  eq.  (2.2.9)  into  (2.1.11)  gives 


♦  .  —V  I  I  -  /(UP-  ■[XF)dt--lAFtl 

°  Cl-  \  « 

g 


+  t/  7^  -  t}  +  e^Q,a)  + 


20' 


(2.2.10) 


Thus  suppressing  secular  terms  in  t  we  get 


2C 


^ - 

(1-  S) 


2tijkC  _ 


(2.2.11) 


g 


Finally,  substitution  of  eq.  (2.2.8)  for  Q,  and  (2.2.11)  for 


(j)2g(5)  i-nto  equation  (2.2.1);  and  introduction  of  the  new  variable 


2u 


\l>  =  —  eA  lead  to  the  following  nonlinear  SchrUdinger  equation  for  the 
g 

wave  envelope  i|)  . 


ac  ,  ~  c  "  - 

C-  *  2C^ 

g 


e  a. 


2c  a 

8 


},p[2,j;= 


a. 


_2 


—  e  ^  ip 
C  a 

g 


(2.2.12) 


where 


gk3  9th‘*(kh)-10th^(Ich)-H9  _  .g^k  g^^g  .  ^  I 

“l  ■  2a)  '  8th3(kh)  ^20)  2sh(2kh)^  gh-(C^)2 


^2sh(2khT  20)^^  ’  2a)tgh-(C^)2Jch2  (kh) 

O 


kh 


gkC 


k2 


(2.2.13) 


k2 


_S^ 


“2  ~  2a)h  ”  4a)8h(2W0cl?(^)  ^  2aj 


^  . 


9th‘*(kh)-10th2  (kli)+-9 
8th3(kh) 


(2.2.14) 


where  a  =  |'p[ 


(2.2.15) 


Equation  (2.2.12)  will  be  rederived  in  section  (2.3)  by  means  of  the 
Whitham  equations. 


2. 3  Third  Order  Equations  Obtained  from  Whitham  Equations 

We  consider  here  the  case  of  shallow  water  and  recover  the 
nonlinear  SchrUdinger  equation  (2.2.1)  by  using  Whitham' s  modulation 
equations. 

(2.3a)  Modulation  equations 

Considering  the  2-D  problem  of  wave-groups  propagating  over  water  of 
slowly  varying  depth,  the  following  five  unknowns  are  usually  chosen  as 
dependent  variables:  the  wave  amplitude  a,  the  wave  frequency  u  , 
the  wave  number  k,  the  average  water  depth  h,  and  the  current  velocity  U. 
To  determine  these  unknowns  we  start  from  Whitham' s  set  of  modulation 
equations,  Whitham  (1974),  p.  556.  Pseudo-phase  consistency  relation: 

+  T-  [g(ti-h)  + — +  0(e**)]=  0  (2.3a.  1) 

2sh(2kh) 

Mass  conservation  equation: 

lOlihl  +  !_  [h  U  +  4  a2  +  0(g4)]  =  0 

3t  3x  n 

Za 


(2. 3a. 2) 
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Wave-action  conservation,  equation: 


3t 


[4^+  OCe-)] 


+  |jt»'  f-  *■  0(e“)] 


=  0 


(2.  3a.  3) 


Consistency  condition: 


3t  3x  ^ 


(2. 3a. 4) 


-  -  % 

where  a  =  [gk  th(kh)3  is  the  Linear  dispersion  relation, 
o'  =  3cr/3t,  and  e  is  a  typical  wave  steepness,  e  =  0(ak). 


Following  Whitham  (1974),  p.  562  and  including  higher  order 
dispersive  terms,  see  Whitham  p.  526,  (which  arise  from  the  quadratic 
part  of  the  Lagrangian  L  =  oj  -a),  the  dispersion  relation  is 
given  by: 

m  =  a  +  kU  +  , —  (h-h)  +  a^  -  —  +0(e'*)  (2.3a.5) 

2ach^(kh)  2a  ^3x2 

here  C  =  (9th‘*(kh)-  lOth^(i^)  +9)/8th3(kh) 


Here,  and  in  what  follows,  we  assume  that  the  small  modulation 
parameter  e  is  of  the  same  order  of  magnitude  as  the  typical  wave 
steepness.  We  also  assume,  as  we  did  in  section  (2.1),  very  mild 
depth  changes,  having  slopes  of  the  order  of  at  most. 


(2.3b)  Induced_Mean  Flow_ 

To  make  the  slow  variation  explicit  cind  to  facilitate  the 
derivation  we  introduce  the  same  multiple  scale  variables,  t  and  5 
as  given  in  section  (2.1)  by  (2.1.5). 


Rewriting  eqs.  (2.3a.  3)  and  (2.3a. 4)  with  the  new  coordinates 
(2.1.5)  and  averaging  them  over  t  gives: 

-  V 

SL  constj  =  B;  oj  =  const^  =  o)  (2.3b.l) 

rr 

Here  we  assume  that  the  behavior  of  the  solution  as  a  function  of  x 
is  the  same  as  that  of  the  linear  boundary  condition  at  x  .  Namely, 

<30 

decaying  for  “  in  the  case  of  wave-packets  and  with  a  constant 

period  in  the  case  of  modulated  wave-trains. 

Again,  bars  Indicate  averaging  over  the  appropriate  domain  in  t 
(finite  for  modulated  wave-trains  and  infinite  for  wave  packets) 
and  B,  uj  are  the  averaged  wave-action  flux  and  the  so-called  carrier 
frequency,  respectively.  Note  that  for  wave-packets  B  *  0. 

The  carrier  wave  number  k,  satisfies  the  dispersion  relation 
(2.1.25). 

Rewriting  eqs.  (2.3a. 1)  and  (2. 3a. 2)  with  the  new  independent 
variables,  yields: 


(2.  3b.  3) 


Naglectlng  the  second  terms  in  the  above  equations  for  the  time 


being,,  we  obtain: 


9  h  G,+C„G, 

^ - + - i — B_£ 

2^(2kh)  2a  gh-(C  )2  gh-(C  )2 

8  8 

akC  ,  8G2+C  G. 

u»-(4r^  +  — - - a_)  - a - .S.  \ 

2a  2sh(2kh)  gh-(C  )2  gh+(C 

o  S 


(2.  3b.  4) 


(2.  3b.  5) 


where  Gj^  and  are  functions  of  C  ,  which  emerged  as  a  result  of  the 
Integration.  Now,  averaging  eqs.  (2.3b. 2)  and  (2. 3b. 3)  and  substituting 
(2. 3b. 4)  and  (2. 3b. 5)  yields: 


[ 

2sh(2kh) 


___S _ 

gb-(Cg)2 


2sh(2kh)  2u> 


h  Gj^+C  G2 

4.gC  - 6 -  const. 

®  gh-(Cg)2 


+ 


(2. 3b. 6) 


h 


gkC^ 


2sh(2kh) 


-)] 


gG  ^-G  G 

^^g  ih-TcV  ”  (2.  3b.  7) 


where,  consistent  with  the  order  of  approximation  k  and  a  have  been 
replaced  by  k,u. 


Following  Stiassnle  and  Peregrine  (1980)  we  assume  zero  averaged 
mass  flow  (.thus,  const^=0)  and  choose  such  a  reference  level  that  h  -  h 
In  deep  water  (which.  In  turn,  sets  const^  =0).  Having  fixed  these 
two  constants  we  solve  eqs.  (2. 3b.  6)  and  (2. 3b.  7)  for  G^^  and  G^  and  then 
return  to  eqs. (2. 3b. 4) ,  (2. 3b. 5)  to  obtain  the  final  results  for  the 


Induced  mean  flow: 
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gkC, 


U(x  o  =  (^  +  --ir S  .  )  .  a-  - a^,  _ 

4  2sh(2kh)^  gh-(C 

g 


gka 

2ajh 


ka^ 


2sh(2kh)  "^  2  gh-(C  "  2sh(2kh) 

g 


(2.3b. 8) 


(2. 3b. 9) 


Consistent  with  our  level  of  approximation  we  have 


^  ,  for  modulated  wave-trains 

Zoic  “ 

^8 


(2. 3b. 10) 


0. 


for  wave-packets 


(2.3c)  Nonlinear  SchrBdlnger  Equation 

To  derive  the  NLS  we  reqrite  eq.  (2. 3a. 3)  as  follows 

2a  3a  o* a^  3k  ^  2q* a  3a  a* a^  3q  ^  3g'  ^  ^ 

o  3t  o  3t  q  3x  q^  3x  q  3x 


Applying  eq.  (2. 3a. 4),  which  gives  mutual  cancellation  of  the 
second  and  fourth  terms  in  the  above  equation,  and  dividing  by 
2a/q  yields 


1?^“' 


3a  ^  3q* 
3x  2  3x 


0 


(2.  3c. 2) 


The  Taylor  aeries  expansion  of 

“  -  (j'l 

q  =  01  +  C  •  (k-k)  +  -r  • 
g  2 


a(k)  in  the  vicinity  of  k  =  k  is 
(k-k)^+0(G2)  (2. 3c.  3) 


{ 


Using  this  series  we  rewrite  eq.  (2. 3c. 2)  as  well  as  the  dispersion 
relation  (2. 3a. 5) : 


9a 

9t 


(u'*a  9  (k  -  k) 

2  ax 


+  a 
2  9x 


(2.  3c. 4) 


0)  =  0)  +  Cg  •  (k  -  k)  +  (k  -  k)^  +  Oja^-  ^  +0^^  (2. 3c. 5) 

where 

and  are  given  in  section  (2.3)  by  eqs.  (2.2.13)  and  (2.2.14) 
respectively. 


Referring  to  eq.  (2.3a. 4)  we  define  a  phase  function  'f  so 


that 


9t  ’ 


(2. 3c. 6) 


Substituting  eqs.  (2. 3c. 6)  into  (2. 3c. 4)  and'(2.3cr5)  we  obtain 


9a 

9t 


+ 


C 

g 


9x  2 


1!?  !§)+  ±  = 

9x  9x''  2  9x 


0 


(2.3c.  7) 


0 


(2. 3c.  8) 


Multiplying  eq.  (2. 3c. 7)  by  the  imaginary  unit  i,  adding  to  it 

iV 

(-a)  times  eq.  (2. 3c. 8),  and  then  multiplying  the  sum  by  e  we  get 

3  C  11  •• 

2*93^  •’1^  +  iC'I't  +  Cg  ip^)  +  txx"  “il'l'P’l'  “  (2.3C.9) 

i(fl 

where  i(i  *  a  e  is  a  complex  wave  envelope.  Alternatively,  using 
the  scaled  coordinates  (2.1.5)  we  obtain  the  N.L.S.  equation 


(2.2.12). 


3.  SOLUTIONS 


In  this  chapter  we  solve  the  third  order  evolution  equations 
for  two  particular  Important  cases.  The  relevant  eqs.  were  derived 
in  chapter  2,  and  are  given  here  again  for  the  sake  of  clarity. 

From  the  Laplace's  equation  (2.1.1)  it  follows  that 


aad 


3 


(3.7) 


In.  section  (3.1)  we  consider  the  case  of  periodical  boundary 
conditions,  and  present  an  approximate  analytical  solution  for  the 
shoaling  of  wave  groups  over  a  very  slowly  varying  topography.  In 
section  (3.2)  we  consider  the  case  of  infinitely  deep  water  and 
calculate  the  induced  mean  flow  accompanying  an  envelope  soliton. 


(3.1)  Shoaling  of  wave-grougs 

Restricting  the  discussion  to  cases  for  which  the  complex 
wave  envelope  ij/(T,5)  is  periodic  in  t  ,  and  assuming  zero  averaged 
(over  t)  mass  flow  in  the  x  direction,  enables  the  decoupling  of 
equations  (3.3)  and  (3.4)  (see  Appendix  A).  For  this  case  is 

governed  by  the  nonlinear  SchrUedinger  equation  (2.2.12): 

_2 

.  3C  _  .  C  "  .  e  a,  -  . 

8  ^  g  8 

“2  e-2-j  I 
=  —  ^  ^ 

8 


with  Oj  and  02  given  by  (2.2.13),  and  (2.2.14)  respectively.  The 


induced  mean  flow  potential  is  given  by 
*^10  ' 


2TtinT 


2<;jh 

-Zninr 


+  b_^(c)  e  Y 


T+  1  a  (z){  b  (5)  e  + 

,  n  n  ^ 

n=l 


+  Q2(C) 


(3.1.1) 


uiiere  Y  Is  the  period  of 


a  =  cosh(2Tren.Cz+h) /C  y) 
n.  g 


(3.1.2) 


b 

n 


_2  2 
e  g  B, 


2..(  ^ 


_  (^)2 


,  .2Trnehv 
cosh(  TT-— ) 


tl 


(3.1.3) 


and  are  the  Fourier  coefficients  of  ([ip|^  ) 


2uinT 


T 

-  2TrinT 


=  t  {C^(Oe  ^  +  C_^(5)  e 

n=l 


(3.1.4) 


The  convergence  of  the  Fourier  series  (3.1.4)  was  assumed  to 
be  independent  of  e  . 

The  potential  (|)2q  which  is  needed  to  calculate  the  mean 
water  level  5  ,  see  eq.  (2.1.27)  is  given  by 


***20 


gkC 

2u)h 


(3.1.5) 


All  technical  details  are  given  in  Appendix  3. 


A  sl-ipler  and  dimensionless  form  of  eq.  (2.2.12)  is  obtained 
by  means  of  the  transformations 


^  = 


2u^C 


^  il)  exp  1  (ji|;p 


/ 


X  a»dx 


(3.1.6) 


wriich  give 


i'Px  '^TT 


(3.1.8) 


-g^C 


(3.1.9) 


The  dimensionless  parameter  y  Is  a  monotonic  Increasing  function  of  kh, 
having  the  values  zero  and  (oiy)^  for  k.h=1.363  and  kh  -*■  “ 
respectively.  The  statement  of  the  mathematical  problem,  given  by 
Eq.  (3.1.8),  is  completed  by  the  following  input  condition  at  X  =  0. 

(i.e.  X  =  X  ;  a  reference  point  in  infinitely  deep  water). 


4i(T,0)  *  1  +  26e  cos(2TrT) 


(3.1.10) 


which  corresponds  to  a  system  composed  of  a  carrier-wave  and  a 
symmetric  "side-band"  disturbance. 

*  00  #il* 


g^-i[(l-2xG/(UY)a)t-a]j^(^2j 


(3.1.11) 


For  constant  depth,  y  =  const,  it  is  well-known  that  Eq.  (3.1.8)  with  T 
in  (0,1),  subject  to  periodic  boundary  conditions  has  the  following 


X  Invariants 


(3.1.12) 


J,  =  /  ('i>*iJ'^-W^)dT 


(3.1.13) 


J.  =  [  (Ul**-  -[il)^12)dT 


(3.1.14) 
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These  Invariants  are  deteimiined  by  the  Input  condition,  Eq.  (3.9)  so 
that 

Jj  =  1  +  2&^;  =  0  (3.1.15) 

J3  =  1  +  (4-P+2cos2o)*2b2  +  6$**  (3.1.16) 

where 

P  -  Sir^/u  (3.1.17) 

For  varying  depth,  p  *  y (X) ,  J ^  and  J2  remain  invariant  and  are  given  by 
(3.1.12),  (3.1.13)  and  (3.1.15),  (3.1.16),  but  is  a  function  of  X 
governed  by  the  equation, 

dJ  2y  1 

■jx"  = -;7  /  (3.1.18) 

(See  Appendix  5) . 


(3.1a)  Three-Waves  systems 

The  solution  of  Eq.  (3.1.8)  can  be  expanded  in  a  Fourier  series 

00 

^(T,X)  =  Z  D  (X)e^^^’’^  (3.1a.l) 

n 

n*— “» 

The  boundary  condition  at  X  “  0,  Eq.  (3.1.10),  gives  D^(0)=1; 

D,(0)  =  D  ,(0)  =  ae^“;  D  (0)  *  D  (0)  =  0  for  n  ^  2. 

1-1  n  -n 

Stiassnie  and  Kroszynsld.  (1982)  truncated  the  above  given  series  and 
considered  only  three  waves  systems ; 

<(;(T,X)  =  E  D  (X)e^’^^^  (3.1a.2) 

n— 1  ^ 

Substituting  Eq.  (3. la.2)  into  Eq.  (3, 1. 8)  yields  the  following  system 
of  ordinary  differential  equations: 


(3. la. 3) 


^  dX  wCcIDoI^+aId^I^)D^+2d2d*3  -  0, 

dD,  p 

i -^  +  y[(2lDoP+3lDi[^-  |)Di  +  dJd*]  =0.  (3.  la. 4) 


Note  that  Eq.  (3.1.13)  yields  D  =  D^. 


For  constant  depth  the  system  of  Eqs.  (3. la. 3),  (3.1a. 4)  has 
exact  solutions  in  terms  of  Jacobian  elliptic  functions  with  periods 
of  order  1  in  X  which  is  summarized  in  Appendix  B;  for  details  see 
Stiassnie  and  Kroszynski  (1982).  These  solutions  depend  on  the 
invariants  and  on  the  parameter  u  ,  which  in  turn  depends  on  the 

water  depth  h  and  on  the  modulation  period  y  •  For  very  mild  depth 
variations,  where  h^  =  o(l),  we  apply  an  asymptotic,  WKB  related 
approach,  assuming  the  local  solution  to  be  that  of  the  constant 
depth  type  and  using  Eq.  (3.1.18)  to  determine  J^.  Jj^  and  y  are 
fixed  by  the  input  conditions  and  J^,  is  given  through  I(P)  by: 
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(3.1a.5) 


(see  Appendix  4) ,  where 

I(P)  =  J3(P)  - 


(3.1a.6) 


The  initial  value  of  I,  at  X  =  0,  where  P  =  P  is  denoted  by  I  and 

o  o 

is  given  by 


I  = 


+  4(l+cos2a)-2P  ] 


(3.  la. 7) 


as  well  as  I(P)  were  assumed  to  be  of  o(l)  thxougkout  the  rather 
lengthy  derivation  of  Eq.  (3.1a.5).  In  all  our  examples  we  choose 
(y  =  Ittuj  correspondiag  to  the.  fastest  growth-rate  of 

the  Benlamln-Felr  instability. 

(3. lb)  Numerical_Verlf ication_of  the_asymgtotic  solution 

In  order  to  appraise  the  relevance  of  the  asymptotic  solution 
given  in  the  previous  section,  we  compare  i^s  results  with  those  of 
a  numerical  solution  of  the  system  of  ordinary  differential  equations 
(3. la. 3),  (3.1a.4). 

Fig.  1  shows  I=I(P)  for  four  initial  values  of  I  =1(2) 

o 

(I^=  -O.OA,  -0.01,  0.04  and  0.1).  The  broken  line  represents  the 
asymptotic  solution  and  was  obtained  by  numerical  integration  of 
Eq.  (3.1a. 5).  The  solid  line  was  obtained,  by  substitution  of 
the  results  obtained  from  a  numerical  solution  of  the  system  of 
O.D.E.  (3. la. 3) , (3. la. 4)  into  the  expression 

I(P)  =  2lDj,p{|Dj,  p+2lD^l^-P+2lD^[2cos[2(arg  D^arg  D^)]}  (3.  lb.  1) 

The  numerical  solution  of  the  system  of  O.D.E  was  obtained 
using  a  trapezoidal  method  and  assuming  the  P(h(X))  =  2+0. 2X. 

Note  that  the  assumption  ll[  «1,  which  is  necessary  for  the 
asymptotic  solution  to  the  valid,  imposes  a  restriction  on  the 
range  of  variation  of  P  (P=2.8  corresponds  to  cd^h/g=4). 


In  Fig.  2  we  show  three  parts  of  the  exterior  group  envelope 

kco.x)!  as  well  as  the  interior  group  envelope  for  the 

input  conditions  a  =  0,  -  0.158  (I  =  0.1). 

o 

Here  again,  solid  lines  represent  the  numerical  solution 
of  Eqs.  (3. la. 3)  (3. la. 4)  with  P  =  2+0. 2X  while  the  broken  lines 
correspond  to  results  obtained  by  the  asymptotic  method,  utilizing 
the  relation; 

[4i(T,X)p  =  {-2I+2(4Jj-P)z-7zMs(2I+2Pz-z2)*^  + 

+  4zcos2irT]^}/8z  (3. lb. 2) 

where  z  is  given  in  Appendix  B. 

The  thre  parts  shown  in  Fig.  2  are  for  P  =2,  2.44,  2.75 
for  the  asymptotic  solution  compared  to  P  in  (2,2.03),  (2.44,2.49), 

(2.75,  2.81)  for  the  numerical  solution  of  the  O.D.E  respectively. 

The  agreement  between  the  two  methods  of  solution,  as  seen  in  both 

the  above  figures  is  rather  encouraging  and  seems  to  indicate  the 

validity  of  our  new  asymptotic  solution  of  the  system  (3. la. 3) (3. la. 4) 

Nevertheless,  one  still  has  to  answer  the  question  if,  and  to  what 

extent,  the  system  (3. la. 3,4)  itself  is  a  reasonable  substitute 

for  the  N.L.S.  equation.  For  constant  depth  Stiassnie  and 

Kroszynskl  (1982)  show  a  good  quantitative  agreement  in  the  length 

of  the  modulation-demodulation  cycle  and  only  a  qualitative 

agreement  for  the  amplitudes.  A  similar  trend  can  be  seen  in  Fig.  3,  which 

compares  to  numerical  solutions,  for  for  the  N.L.S.  (3.1.8)  -  dotted 

line,  and  the  other  for  the  system  o'"  O.D.E.  (3.  la, 3, 4)  -  solid  line. 


The  Input  data  in.  Fig.  3  Ls'^;  =  0,  &  =  0.1  (1^  =  0.04)  and 

the  variation.  P  =  2  +  0.2X  is  assumed;  Both  the  exterior  and  interior 
group  envelopes  are  drawn. 

IJe  believe  that  our  much- simplified  asymptotic  solution  is  not 
over-simplified,  and  is  able  to  produce  quite  a  few  results  of 
qualitative,  and  maybe  even  serai-quantitative  relevance,  which 
enable  us  feome  new  physical  insight. 

(3.1c)  On  P  and  I 
- o 

One  fundamental  property  of  the  asymptotic  solution  is  that 

it  depends  on  P  and  I  solely.  Given  the  input  data  a,  B  (and 

o 

P  =2),  I  is  determined  by  Eq.  (3. la. 7).  Then,  integrating  Eq. 
o  o 

(3.1a. 5)  from  P^  to  P  the  parameter  I(P)  is  found,  and  the  solution 
given  by  Stiassnie  and  Kroszynski  (1982),  is  locally  applied. 

Fig.  4  gives  the  relation  between  P  and  nondimens ional 

local  water  depth  k  h,  (where  k  =  m^/g  is  the  wave-number  in 
^  00  '  00 

infinitely  deep  water) ,  for  y  =  2x0)  .  P  is  a  monotonic  decreasing 

function  having  the  values  infinity  at  k^h  =  1.195  (kh=1.363) 

and  2  for  k  h  The  input  value  1  ,  (for  P  =2)  dependence 

oo  o  o 

on  a  and  B  is  shown  in  Fig.  5.  Note  that  different  combinations 
of  a  and  B  give  the  same  I^,  and  thus  basicly,  the  same  solution 


(3. Id)  Group  enveloges 


The  free-surface  of  an  (unstable)  shoaling  wave-train,  displays 
three  distinct  length  scales:  -  the  wave  length;  -  the  modulation 

or  group  length;  and  -  the  modulation-demodulation,  group-envelope, 
or  maybe  best  'supergroup*  length.  These  three  lengths  are  given  by 


Z^  =  iKfk 


(3.1d.l) 


=■  2ire  foi. 
^  g 


Z 
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-2  3 

4e  P(C  )^ 

oj2aj"/p(4-P) 


In 


2P^(4-Py^ 

(2P-1)I 


see  Fig.  7  a  . 


(3.  Id. 2) 


(3.  Id. 3) 


It  can  easily  be  seen  that  in  the  range  of  depths  where  the  asymptoti 
solution  applies,  k^h  >  4,  and  Z^  remain  almost  constant.  On  the 
other  hand,  Z which  depends  on  P,  exhibits  quite  a  remarkable  variation 
as  shown  in  Fig.  6  . 


In  Fig,  j  ,  we  show  the  variation  of  <2.  ^  as  a  function  of  the 

depth  k^h  for  four  different  input  data  I^=-0.04,  -0.01,  0.04,  0.1.  For 

I  <0,  Z.~  decreases  with  decreasing  depth,  but  for  I  >01, 
o  i  O-i 

increases  with  decreasing  depth  up  to  a  "critical  depth"  (corresponding 
to  I  =  0)  and  from  there  on  starts  to  decrease. 


Fig.  7  shows  the  group  envelopes  (dashed  line)  and  wave 
envelope  (solid  lines)  at  a  fixed  instant  for  e  =  0.2,  at  the  following 


four  Locations:  (a)  -  infinitely  deep  water,  =  2,  =  0.1; 

(b)  =  11.2,  P  *  2.2,  I  =  0.052;  (c)  V.Ji  =  5.7,  P  =  2.45,  I  =  0.01; 

(d)  k.„h  =  4.2,  P  =  2.75,  I  =  -0.028. 

In  Fig.  7  a  we  have  added  a  portion  of  the  wavy-surface  (thi’n 
solid  line)  as  well  as  the  lengths  and  Note  that  the  super¬ 

groups  (namely:  the  exterior  and  interior  group  envelopes)  are  fixed 
in  space,  while  the  wave  envelope  moves  'v^th  group-velocity  and  the 
waves  themselves  with  the  phase  velocity.  Similar  sketches  to  Fig.  7 
were  obtained  for  the  other  cases  given  in  Fig.  6 

In  order  to  complete  the  picture  for  shallower  water  depth  we 
present  in  Fig.  8  the  group  envelopes  at  five  locations: 

P  =  (2,2.05),  b:(2.32,  2.63);  c:(17.7,  22.8),  d:(-22.4,  -17.5), 

e: (-4. 1,-3. 9) ,  as  obtained  from  a  numerical  solution  of  the  system 
(4.3a,b)  for  the  same  input  data  as  in  Fig.  7  assuming  w  »  39.47-lOX. 

The  results  in  Fig.  8  indicate  that  continues  to  shorten 

and  that  the  intensity  of  modulation  decreases. 

(3.1e)  The  mean_flow  field 

We  express  the  mean  flow  U  =  V  =  a<|)^/3z*3'i' /9x 

through  the  stream  function 
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where  Z  =  irez/C  y  .  H  »  ireh/C  Y  and  ’i’  =  e  are  dimensionless 
S  o 

quantities.  The  constant  in  Eq.  (3.le.l)  is  chosen  so  that  =  0 
at  the  bottom.  The  mean  free  surface  is  given  by  Eqs,  (2.1.16) 
and  (2.1.27). 


The  stream- function  ’i'(T,Z)  as  well  as  the  mean  free-surface 
for  cases  a,  b  and  d  of  Fig.  7  are  presented  in  Figures  9  ,  10.  .  and 

11  respectively.  These  figures  demonstrate  the  rather  complicated 
structure  of  the  wave  Induced  mean  flow  field. 


Some  of  the  main  features  are:  (i)  the  mean  current,  which 
is  shown  in  part  (b)  of  the  figures,  as  well  as  the  mean  free  surface, 
in  part  (c)  exhibit  a  somewhat  cellular  structure  influenced  by  the 
wave  envelope  variations;  (ii)  a  dominant  adverse  current  appears 
underneath  the  high  waves  and  a  much  weaker,  positive  current  (in  the 
wave  propagation  direction) under  the  low  waves,  for  the  shallower  cases 
the  positive  currents  almost  disappears;  (ili)  the  magnitude  of  the 


maximum  adverse  currents  at  the  free  surface  is  almost  the  same  for 


all,  three  depths  (k.^K  ■=  11.2  and  4.2);  (tv)  one  can  notice  the 

tendency  of  the  flow  fields  to  become  more  uniform.  In  their  lower  parts 
and  on  the  sides  of  the  supergroups  (where  the  modulation  amplitudes 
get  much  smaller) ;  (v)  there  Is  a  set-down  In  the  mean  free  surface 
accompanying  the  peaks  of  the  wave  envelope  and  a  smaller  set-up 
accompanying  their  troughs. 

(3.2)  The  -induced  flow  accompanying  an  envelope  sollton 

We  consider  here  the  case  of  Infinitely  deep  water.  In  this 

case  we  assume  that  the  induced  mean  flow  4  is  of  order 

o 

(<^10  “  0(e);  (|)2q  =  0(e)).  For  Infinitely  deep  water,  system  (3.1) 
to  (3.4)  becomes: 


ho  ho  '  °  -«  <  z  <  0 

ZZ  g  XT 

(3.2. 

<t'io  “  ° 

z-»— «  z 

(3.2. 

<l>10  =  k([4>12)  ,  Z  =  0 

z  ^ 

(3.2. 

1  ~  _2.._ 

”  g  ’*'tt  '  ^  '(»  =  0,  z  *  0 

(3.2. 

Introducing  the  dimensionless  variables 


Note  that  X,  T,  Z  aad  'I'  ia  the  present  section  are  slighly  different 


from  their  previous  definition. 

to  (3.2,4)  becomes 

=  0 

0 

Z  =  0 
tjf 

0,  z  =  0 


System  (3.2.1) 


A  A 

ZZ  ^  ^TX 


lim  « 


(3.2.6) 


(3.2,7) 


(3.2.8) 


(3.2.9) 


One  well-known  solution  of  equation  (3.2.9)  which  decays  at  T  -  i* 
(x  =  ±«  )  is  the  envelope  soliton,  given  by: 


i|;(X,T) 


=  sech(/2  T)  e^*^^ 


(3.2.10) 


Substituting  (3.2.10)  in  (3.2.8)  yields 


$Z  =  tsech^(.^  T)}  ,  Z  =  0  (3.2.11) 

Equations  (3.2.6),  (3.2.7)  and  (3.2.11)  define  a  Neumann  problem 

in  the  lower  half  plane,  which  was  solved  utilizing  the  Fourier  transform 

method,  with  the  following  result. 

^  00 

=  T  E  (Z-a  )/[X^+(Z-a  (3.2.12) 

n  n 

n=o 


where 


a  =s  —  (!jfn) 


(3.2.13) 


-41- 


All  technical  details  are  given  In  appendix  6.  Let  U  and  V  be  the 
Induced  mean  flow  velocity  components  In  the  directions  x  and  z  respectively. 
From  (3.2.12)  and  (3.2.5)  It  follovs 


u  =-T-  =  ^  (Z-a  )  [CZ-a  )2  -3X2]  /[T2+(Z-a  )2]3 
€  01  ti  n.  n 

n*0 


(3.2.14) 


V  =  pjij  =  T  Z  [T2-3(Z-a^)2]/[T2+-(Z-a^)2]3 
n=0 


(3.2.15) 


The  dimensionless  stream  function  is  given  by 


2  «  T2-(Z-a  ,2 

’»'(x,z)  -■V'F  =  %  +  %  Z  - - - 

^  “  n=0  [T2+(Z-a  )2]2 

n 


(3.2.16) 


Note  that  1  ®n~  »  1  and  that  the  first  tern  on  the  r.h. s.  of 

n=0 

eq.  (3.2.16)  was  chosen  to  render  'i'(0,0)  ■=  0.  The  streamlines  of  the  induced 
mean  flow  field  are  shown  in  Fig.  12,  the  total  flux,  per  unit  width, 
involved  in  this  flow  is  given  by 


q  »  ^(-.O)  -  ^(O.O)  =  %  €2  ^ 


(3.2.17) 


This  value  is  equal  to  the  Stokesian  mass  transport  at  the  peak 
of  the  wave  packet,  as  it  should  be. 


•  .*  •.•] 


4.  CONCLUDING  REMARKS 


The  main,  result  of  the  present  study  is  the  development 
of  an  analytical  solution  which  is  able  to  provide  detailed 
information  about  the  physical  quantities  involved  in  the 
shoaling  of  wave  groups. 

The  quantities  Include:  the  variations  of  the  mean  free 
surface  (also  called  set-up  and  set-down);  the  induced 
mean  flow;  and  the  supergroups  (group  of  groups)  which  to  our 
knowledge  have  not  been  discussed  in  the  past. 

The  improvement  in  our  understanding  of  the  above  mentioned 
quantities  is  of  practical  importance  since  it  is^  believed  that 

i 

they  are  related  with  such  phenomena  as  surf-beats,  longshore  : 
cellular  structure  and  Harbour  resonance. 


conditions 


If  ill  and  <|i^Q  are  periodical  of  period  Yf 
to  (3.7)  is  given  by: 


the  solution  of  (3.  1) 


-i^  /  dx 

i<C,T)  -  ®  0^(0  + 

O' 

o 

"  2xlnT  -Zirint 

+  ^  {D  (5)  e  ^  +  D  (5)  e  ^  } 

,  n  "n 


(A.l) 


where  a_  is  given  by  (2.2.14). 


2TiinT 
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<&10(C,T,z)  =  Q^(Ot  +  Q2(0  +  ^  a^(z){  b^(Oe  + 

n=l 

-2irinT 
+  b_^(5)e  ^  } 


(A.  2) 


Substituting  (A. 2)  in  (3.1),  and  using  the  b.c.  (3.2)  yields 


en'2Tr(n+h) 


a  =  cosh[ - ;; - 

n  C  Y 


"3  ,  n-  1,2,... 


(A.  3) 

(same  as  (3.1.2)) 


From  (A. 3)  and  the  b.c.  (3,3)  we  obtain 


_ _ 


_  e  ^g^^2 


(A.  4) 


eV 


C  Y  (same  as  (3.1.3)) 


where  C  (f)  are  the  Fourier  coefficients  of 

n  ^  T 


2  . 


2iTinT 


-2TvinT 


t  (Same  as  (3.1.4)) 


Substituting  (A. 3)  and  (A. 4)  in  (A. 2)  and  differentiating  twice 


with  respect  to  f  yields: 
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Thus  from  (A.  5)  and  (A.  6)  follows  the  relation 
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TT 
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(A.  8) 
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(A.  9) 


r2^ 


It  can  be  easily  seen  that  rQ“rj  <  0(e  )  for  all  n.  On  the  other 


hand  C  (O  are  Fourier  coefficients.  Thus  there  exists  a  N,  such  that 
n 


2TinT 


-2irinT 


Z  (0(0  e  '  -hC  „(5)e 

XT 

n=N 


j-n  Y  J  ^  g.  3 


and  then 


(A.  lo: 


Zirinr  -27rinT 


Z  (r  -r.)  {Ce^  -t-Ce  ^1<  0(e) 


For  n.  <  N: 


r  -r,  <  Q(e"b 

n  1 


From  (A.  7),  (A.  10)  aad  (A.  11)  we  finally  obtain  the  relation 


'10 


TT 


=  +  Q(e) 


z=0 


And  using  the  relation 


.-2  2 

^  8  02  -1 
r. - - +  0(e  n 


^  2(0(1  -  1^) 


which  can  be  easily  verified,  we  obtain 


10 


TT 


=  rl(p|2+  0(e) 


z«0 


where 


.-2„2 


e  g 


r  = 


2(0  (1  -  fy) 


(A 


(A, 


(A, 


(A, 


(A 


Equation  (A.  14)  is  identical  to  (2.2.4)  (A  •  (p  ),  which 


was  obtained  for  the  case  of  shallow  water.  Integrating  (A. 14)  with 
respect  to  t  yields 


10 


=  rjl^P  +  Q(5)  +  0(e) 


(A 


z=0 


From  here  we  proceed  in  a  similar  way  as  in  section  (2.2)  and 
obtain  that  Q(0  Is  given  by  (2.2.8),  by  (2.2.11),  and  Q,  (0  =  i}> 


Finally  aubstltutlag  ,  Q(S)  and  <{>20  in  (3.4)  we  obtain 
that  the  wave  envelope  i|)(C,?)  is  governed  by  the  nonlinear  SchrHedinger 
equation  (2.2.12). 

^  _  C"  _  e~^a,  .  ,  e~^o  —  . 

^^ip  +  iij.^  +  — 0-^  kP'P  =-0 — \^\H 

g  g  g  g 


where  and  02  given  by  (2.2.13)  and  (2.2.14)  respectively. 


For  deep  water  (eh  »  YC^)  ,  eqs.  (A..  8)  and  (A..  9)  become 


(A.  17 


(A.  18 


respectively,  and  one  can  easily  verify  that  (}),„  =  0(e). 


r,  = 


e" 

‘^Yg 


and 


r  =  - 

n 


-1  y 
ojyg 


Input  data: 


-i  (4J,-P){I-[l-jJ^]  ’'-) 

J,  -  P  [l-(l+%'^] 

^  P 

c  =  max(cj^,C2) ;  ~  mln(Cj^,C2) 

e  =  2P;  b  =  y(4-P) 

^2  _  1  e-b  .  c-d 
e  b 

/Teb  •  X 

y  =  ._ 


(B.l) 

(B.2) 

(B.3a,b) 

(B.3c,d) 

(B.4) 

(B.5) 


cd(y,r)  is  a  Jacobian  elliptic  function  of  the  argument  y  with  modulus  r. 


ebd-cd^) 

z  =  — r - 77 — ~ 

e-b*  cd 


(B.6) 


IdJ  =  /I72  ;  1d^I=  /.Tj-z 


(B. 7a, b) 


I  -1 . 5z^  -*-(P-2J  1 )  •  z 
2z(J^-2) 


cos[2(arg  D,-arg  D  )]  = 


(B. 7c, d) 


Appendix  C  -  Derivation  of  the  Equation  for  J^CX) 


Let  us  suppose  that  u  Is  a  given  function  of  X.  Differentiating 
equation  (3.1.8)  with  respect  to  T  and  multiplying  by  we  obtain 


Subtracting  from  (C. 1)  its  complex  conjugate  and  dividing  the  result 
by  p(X)  yields 


1  1  3 

Hex)  *  ~ax  u(x)  "§!  =  o 


(C.2) 


Multiplying  equation  (3.1.8)  by  t(<*,  taking  the  complex  conjugate 
and  substratting  it  from  the  resulting  expression,  we  obtain 


1  al  =  0 


3T 


(C.3) 


Multiplying  (C.3)  by  yields 


1*1^  (V*-*!*)  -» 


(C.4) 


Subtracting  (C.4)  from  (C.2)  gives 


^  (C.5) 


y(x)  ax  2  ax  ai 


which  can  be  written  as 


I'P. 


.  d  r  ■  i.-  ,  I  ,  lUi  .  ■  i-  .  9  r  ^T'^TT'  . 

‘  '  Too  ■  ^1*1  ’ 


(C.6) 


Integrating  (C . 6)  with  respect  to  T  from  0  to  1 ,  and  using  the 


fact  that  is  periodic,  with  period  1  we  obtain: 


3X  ■’  U(X)  2 
o 


J  I  r.j.  I 
O 


(C.7) 


which  can  be  written  as: 


^JjCX)  Ifil^dT 

O 


(C.8) 


where  J^CX)  =  /  ■  |<J/^l2)dT. 

Substituting  the  approximate  solution  (3. la. 2)  into  the  right 
hand  side  of  (C.8)  we  obtain  the  following  approximate  equation  ^or 


-Tx-  - 


(C.9) 


where  z(X)  =  2|Dj^P. 


Integrating  eq. (C.9)  yields 


J-(X)  =  J-(0)  +  8Tr2/*  ^  z(X)dX  =  J  (0)  + 


(C. 10) 


and  using  (3.1.17)  one  finally  obtains 

P 

J,(P(X))  =  Jo(P„)  -  /  z(P)dP 


(C.ll) 


Equation  (C.ll)  is  an  integral  equation  for  J2(P).  If  we  suppose  that 

2 

J^(P(X))  varies  in  such  a  way  that  then  z(P)  may  be  replaced  by 

the  solution  for  constant  depth  (B.6)  so  that 


J  (P)  =  J  (P  )  -  2P14^)_(1_-Cd:_) 

s'-  ^  s'-  7P-(4-P)cd^ 


where  cd  is  the  Jacobian  ellyptic  function  with  argument 


and  modulus  r  given  by 


2  _  ,  (e-b)(c-d)  _  ,  ,  (4P-2)  2,2 

1  eb  ^  ^  p^(4-P)^  ^“^S  0(23-^^ 


e,b,c,  and  d  are  given  in  Appendix  B. 


(C.12) 


(C.IS) 


(C.14) 


After  some  algebra,  4hd  using  the  relations 


J  cn 
cd  =  — 
dn 


2  2 
cn  =  l-sn 


(C.15) 


.2  ,22 

dn  =  1-r  sn 

,2  ,  2 

r  =  1-r 

where  cn,  dn  and  sn  are  Jacobian  ellyptic  functions,  equation  (C.12) 


is  written  in  the  form 

J  (P)  =  JoCP^)  -  / 

3  J  o  p'  e-D 

where  ^ 

_2  _  ,  (c-d) 


1-q^sn^ 


(C.16) 


(C.17) 


The  integrand  in  (C.16)is  an  oscillatory  function  of  the  argument  y, 
and  if  we  again  apply  the  very  mild  topography  assumption,  it  is  justified 
to  replace  it  by  its  average  as  follows 


l-q'^sn‘^(y)  e-b 


l-q-)r‘ 


where  4k(P)  is  the  period  of  the  sn  Jacobian  ellyptic  function, 
is  the  Jacobian  Zeta  function  and 

P  =  sin  ^(-J) 

In  order  to  obtain  (C,18),  equation  (414.02)  from  Byrd  and 
Friedman  (1971)  was  used. 

Substituting  (C.18),  into  (C. 16)  one  obtains,  after  some  algebra 


P 

J3(P)  =  J3(p,r)dP  (C.20 

P 

o 


or 

ap-  =  Zj(p,r) 

(C.21 

where 

sln^p  “  1  -  ^  +  O(S^) 

(C.22 

5) 

eb 

(C.23 

6  (c  d)  p(4-p)  ^■^3  ^ 

(C.24: 

Finally,  using  the  fact  that  sin^p  = 

1-H)(d),  and  r^=l-H)(6),  we 

replace  Zj(p,r)  in  (C.21)  by  a  simpler  expression  as  follows. 

The  function  Zj(p,r)  is  defined  by 

Zj(p»r)  =  E(p,r)  -  I  •F(p,r)  (C.25) 


Here  E,  K,  E(p,r)  and  F(p,r)  are  the  Ellyptic  Integrals.  (See  Byrd 
and  Friedman  (1971)). 


The  above  equation  Is  written,  using  CB. 3)  as 


Finally,  Introducing 
I(P)  =  J3(P)  - 


(C.34) 


(C.35) 


(C.36) 


We  consider  the  problem 


(3.2.6) 


lim  $  =  0 


(3.2.7) 


z- 


*z“*^|t  (v'I  T)}, 


z  =*  0 


(3.2.11) 


Equation  (3.2.6)  with  the  boundary  conditions  (3.2.7)  and  (3.2.11) 
constitute  a  Ne\imann  problem,  which  Is  solved  £or  the  Induced  flow. 


Substituting  the  following  Fourier  Transform  of  the  Induced  flow 


potential  function 


f(X,Z)  =  — —  /  e^^''^(X,T,Z)  dT 
/2it  -« 


(D.l) 


Into  (3.2.6),  (3.2.7)  and  (3.2.11)  yields 


£22(1,  Z)  -X'^f(X,Z)  -  0 


-OK  <  Z  <  0 


(D.2) 


1 


£2^.0)  ^  . 

2/2ir  -0 


—Csech^/ZT]  dT,  Z=0 


(D.3) 


11m  f„(A,Z)  =  0 


(D.4) 


f,(X,Z) 


The  solution  of  (D.2),  (D.3),  (D.4)  is  given  by 

1  i^i==r”ixt9 


-  e  /  e^  ^  —  [sech^ii^T  1  dx 


(D.5) 


Then,  by  means  of  Che  Inverse  Fourier  Transform,  we  obtain 
®-(X,T,Z)  e^"^  |-  Lsech^i'^r]  dtdX  = 


^Tl-  (8ech^'^T)[  dT 


The  Integral  In  the  square  brackets  can  be  evaluated,  as  a  sum  of  two 


Integrals,  from  -«  to  0,  and  from  0  to  “  ,  yielding  the  result 
1  ■*“ 

-  -  2"  L  •f^Csech^.^T]  dt  <“•’ 


Finally,  integration  with  respect  to  Z  gives 
~  -1-  s 

4(5^T,Z)  »  -4ir  /  ln[(T-T)^4-Z^]‘|:^C8ech2/2T]dT+  Q(X,T) 


From  equation  (3.2.6)  it  follows  that 


Q(X,T)  -  a(X)T  +  W(X) 


and  imposing  the  condition 


11m  =  0 

T-+±oo  ^ 


yields  that  Q  is  only  function  of  X. 


Q  »  W(X) 


On  the  other  hand,  replacing  x  (D.8)  by  the  new  variable 
x-T,  and  integrating  twice  by  parts,  yields  the  result 


The  path  T  is  taken  between  two  poles,  as  shown  in  the  figure 


-5 


intBgral  dXong  th.B  half  circla  t6nds  to  zaro  whan  R  , 
and  after  some  algebra  we  obtain  the  result: 

—  -f  Z  ®®  "■  z^ 

f(X,T,Z)  =i|-  {8ech^[/2(T-lZ)]  +-i  .  S  -J -  +  W(X)  (I 

^  i=0  (z2+z2)2 


where  Z.  is  given  by  (D. 13) 
From  the  identity 


sech^C  =  -  Z 


(See  eq,  (3.64)  in  Carrier  et  al.  (1966)),  we  obtain  the  final  result 
^  e^TlZ-  -  (%+k)} 

♦  (X.T.Z)  S  - —  ■  -  +  W(X)  (D.l 

k-0  Ct^  +  (Z--  C'5+k)}2]2 
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1  I*I(P), - Asymptotic  solution,  - Numerical  solution  of 

Eqs.  (3.1a.3),  (3.1a.4). 

2  Exterior-group-envelope  |)|/(o,X)j  and  Interior-group  envelope 
l'l'(%»X)l  for  the  input  conditions  a  =  0,  B  =  0.158,  (1^=0. 1) 

-  Asymptotic  solution,  -  Numerical  solution  of  Eqs. 

(3. la. 3),  (3.1a.4) 

3  Exterior-group-envelope  |)/)(0,X)|  and  Interior-group-envelope 
l'K*SjX)|  for  the  input  conditions  a  =  0,  B  =  0.1,  (1^=0. 04), 

- Numerical  solution  of  O.D.E.  (3. la. 3) , (3. la. 4)  -.-.  Numerical 

solution  of  N.L.S.  (3.1.8) 

4  P-P(K  h)  for  P  =■  2. 

00  O 

5  I  -  I  (a, 6)  for  P  -2  (on  B^  =  0;  I  =0). 

00*  o  o 

6  The  'supergroup*  length  as  a  function  of  depth  h. 

7  The  group  envelope  ( - )  and  wave  envelope  ( - )  at 

t  =  constant  for  P  “2,  I  =0.1,  e  =  0.2  and  K  h  =  (a)“>  , 
o  o 

(b)  11.2,  ^c)  5.7  and  (d)  4.2. 

8  The  group-envelopes,  for  Pq“2,  1^=0. 1,  at  K^h  =  (a)®,  (b)7. 11-4.8, 

(c) 1.32-1.29,  (d)l,12-l.ll,  (e)0.96-0.95. 

9  The  flow  field  for  1^=0. 1,  e=  0.2  at  K^h  =  ®.  (a)  group- 

envelopes  and  wave  envelope  (b)  mean  flow  stream-lines,  (c)  mean 
f ree-surfaces. 


10 

As  in 

9 

for  K  h  =  11.2 

oo 

11 

As  in 

9 

for  K  h  =  4.2 

00 

12 

The  wave- 

induced  mean  flow  field  under  a  deep  water  envelope 

soliton 


Exterior  -  group-envelope | i|i (o,x)  ]  and  Interior-group  envelope [ >!j (Jj.x)  |  for  the  Input  conditions 
o  »  0,  6  *0.158,  (I  *0.1) - Asymptotic  solution,  - Numerical  solution  of  Eqs.  (3.  la.3  ) 


*0.04),  - Numerical  solution  of  O.D.E.  (3. la. 3),  (3. la. 4),  - -  Numerical 


The  group  envelope  ( - )  and  wave  envelope  ( - )  at  t-conatant  for 

V2»  e-  0.2  at  k,h-(a)«,  (b)  11.2,  (c)  5.7  and(d)  A. 2. 


6-0-96 


Exterior-group  envelope  (a) 


- 1 - i ^ - i - ^ - 1 - 1 - 1 - 

T=Aiie”^X 

Figure  9 

The  flow  field  for  ^^•2,  e  *0.2  at  ■  *  (a)  group  envelopes 

and  wave  envelope  (b)  mean  flow  at ream- lines,  (c)  mean  free  surface. 


t-lnduced  mean  flow  field,  under  a  deep  water  envelope  sollton 


